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Abstract 

We derive Thouless- Anderson-Palmer (TAP) equations for quantum disordered 
systems. We apply them to the study of the paramagnetic and glassy phases in 
the quantum version of the spherical p spin-glass model. We generalize several use- 
ful quantities (complexity, threshold level, etc.) and various ideas (configurational 
entropy crisis, etc), that have been developed within the classical TAP approach, 
to quantum systems. The analysis of the quantum TAP equations allows us to 
show that the phase diagram (temperature-quantum parameter) of the p spin-glass 
model should be generic. In particular, we argue that a crossover from a second 
order thermodynamic transition close to the classical critical point to a first order 
thermodynamic transition close to the quantum critical point is to be expected in 
a large class of systems. 
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1 Introduction 



Glassy systems of extremely diverse types exist in nature. They all share several common 
features like a very slow, non-equilibrium dynamics. The development of a full theoretical 
description of the glassy phase is one of the most important challenges in condensed 
matter physics. A variety of techniques, that range from scaling arguments to mean- 
field approaches have been, and are still used, with the aim of attempting a satisfactory 
description of the glassy properties. 

One of these techniques is due to Thouless, Anderson and Palmer (TAP) |^, who 
introduced an approach to classical disordered systems based on the study of a free-energy 
landscape. The key object is the Legendre transform of the free-energy F{(3) = — \nZ/(3 
with respect to a number of order parameters that are sufficient to describe the transition 
and the different phases in the system. This function behaves as an effective potential 
whose minima represent different possible phases. In a classical fully-connected Ising 
model only one order parameter is needed, the global magnetization m = J^ii^d/^- 
The two possible minima of F{(3,m) correspond to the two possible states of positive 
and negative magnetization, m = ±mo(T). Focusing on the Sherrington-Kirkpatrick 
(SK) mean-field model for spin-glasses, TAP showed that all the local magnetizations 
fTT-i = (si), i = 1, . . . . , N, have to be included in order to derive the relevant free-energy 
landscape. The extremization condition of the TAP free-energy on the mjS leads to the 
TAP equations. It was soon after realized by Bray and Moore [0, 0] that the number 
of solutions to the TAP equations for the SK model is exponential in the number of 
spins in the system, for temperatures below the spin-glass transition [Q. A very useful 
alternative derivation of the TAP equations was given by Plefka ^ who showed that 
these equations can also be obtained from a power expansion of the Gibbs potential up 
to second order in the exchange couplings. The advantage of this derivation is twofold: 
it allows to show convergence of the power expansion for all temperatures and it is easily 
applicable to other mean-field glassy models. Moreover, Georges and Yedidia showed 
that the high temperature series, at fixed order parameter, of the free-energy can be 
used to derive TAP-like equations, and its corrections, for models in finite dimensions or, 
equivalently, with finite range interactions. The connection between the TAP approach 
and the more standard analysis of the partition function of a disordered model has been 
exhibited by De Dominicis and Young 0] who showed that, for the SK model, one recovers 
the equilibrium results of the replica or the cavity method via weighted Boltzmann 
averages over solutions of the TAP equations. More recently, the TAP approach has 
been applied to other classical disordered models. In particular, two models that we shall 



discuss in the following, the spherical and Ising p spin-glass models [^, |T0|, |TT], [T^, |T^, [14 1 
and the Ghatak-Sherrington (GS) model [^, |T^ have been analyzed with this method 

EH SID. 

Glassy systems, and in particular disordered ones, are characterized by having a very 
slow dynamics with non equilibrium effects at low temperatures [pO|, Mean-field 



models, like the the spherical p spin-glass model |^ or the SK spin-glass [^, capture 
this phenomenology. The dynamic solution for the evolution starting from random initial 
conditions, that represent a quench from high temperatures analytically, is intimately 
connected to the structure and organization of TAP solutions. One of the most striking 
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results of the dynamic analysis of p spin-glass-like models is that the energy density (and 
other one time-quantities) converges asymptotically to the energy density of high lying 
solutions of the TAP equations. This level has been called threshold. The energy-density 
in equilibrium is different. This and other related results suggest that an interpretation of 
the dynamics in terms of a motion in a TAP free-energy landscape can be given . The 
generalization of the TAP approach to dynamics that has been developed in p4| allows 
one to make this statement precise: the evolution is determined by a gradient- descent 
in the TAP free-energy landscape with the most important addition of non-Markovian 
terms. 

Usually, glasses can be analyzed with a fully classical approach since their transition 
temperatures are rather high. Nevertheless, in many cases of great interest the critical 
temperature can be lowered by tuning another external parameter and quantum fluctu- 
ations become very important. This is the case for the insulating magnetic compound 
LiHo^Yi_j^^F4, that is an experimental realization of a quantum spin-glass, and presently 
receives much attention [^. Other examples where glassy properties in the presence 
of quantum fluctuations have been observed are mixed hydrogen bonded ferro-antiferro 
electric crystals [^, interacting electron systems [^, cuprates like La2-xSrxCu04 [p8| . 



amorphous insulators p9[, etc. 

The quantum fluctuations in LiHoa,Yi_3;F4 can be controlled by tuning the strength of 
an external field that is transverse to the preferred direction of the randomly located mag- 
netic impurities. After a series of experiments presented in the authors' conclusions 
are: (1) The samples undergo a paramagnetic to spin-glass transition in the (T, P) plane, 
where P oc and Ht is the strength of the transverse field. (2) The transition is of second 
order (in the thermodynamic sense) close to the classical critical point (T = Tc, P = 0) 
but crosses over to first order close to the quantum critical point (T = 0,P = Pc). (3) 
The system undergoes out of equilibrium dynamics in the glassy phase as demonstrated 
by the fact that the dynamics strongly depends on the preparation of the sample for all 
subsequent times explored experimentally. 

The theoretical study of quantum spin-glasses started with Bray and Moore's anal- 
ysis of the equilibrium properties of the fully connected Heisenberg model [^. In 
this article. Bray and Moore introduced a path-integral representation in imaginary 
time of the partition function that they analyzed with the replica trick. Many arti- 
cles on the equilibrium of this, and related, mean-field models have been published since 



3T| p2| , p3| , |34| , P5| , |36| , |37| , |38| , |39[| . The static properties of low dimensional models have 



been studied and it has been shown that, in finite dimensions, Griffiths-McCoy singu- 
larities are very important close to the quantum critical point |40]. In all these models, 
the transition from the paramagnetic to the spin-glass phase has been reported to be of 
second order throughout. 

In most classical disordered models studied so far the transition from the disordered to 
the ordered phase is of second order in the thermodynamic sense. In the exact solution of 
the SK model, the spin-glass order parameter q{x) is continuous at the transition which is 
of second order in the thermodynamic sense In other classical glassy models like the 
Potts glass 1^ or the spherical |4l| and Ising |^2|, ^ p spin-glasses, the order parameter 
jumps at the transition which, however, is still of second-order in the thermodynamic sense 
since there is neither a jump in the susceptibility nor a latent heat. A classical model 
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that exhibits a first-order transition is the anisotropic p spin-glass, p > 3, in which the 
spins take integer values between —S and 5* and there is an extra term in the Hamiltonian 
—D^isj, proportional to a coupling constant that controls the crystalline tendency. 
In this case, a crossover from a second-order to a first-order thermodynamic transition 
in the plane {T/J,D/J) has been exhibited in the exact solution [Q. The classical 
Ghatak-Sherrington (GS) model [0 is another candidate to exhibit a second to first 



order crossover in the thermodynamic transition. It is the anisotropic extension of the 
SK model, or the p = 2 limit of the previous model. In this case, a crossover from a 
second order to a first order transition in the plane (T/J, D/J) has been exhibited in an 
approximate solution (one step replica symmetry breaking) |T^, The exact solution 
has not been derived yet and it is then not well established if it has a true first-order 
thermodynamic transition. 

In quantum problems, first order transitions have been reported in three models. 
The first one is the so-called "fermionic Ising spin-glass" analyzed by Oppermann and 



collaborators [Q. This model, however, is thermodynamically equivalent to the classical 



GS model discussed above [46|. The other two models are very similar indeed and they 



are different ways of extending the classical spherical p spin- glass model |^ to include 
quantum fluctuations. In one case, the continuous spins are generalized to M component 
vectors and a global spherical constraint as well as commutation relations are imposed 
36] . The other one uses the fact that the spherical p spin-glass model can be interpreted 



as a particle moving in an infinite dimensional hyper-sphere with a random potential. 
Quantization is then done by imposing commutation relations between coordinates and 
momenta fS^, The latter can also be interpreted as an extension of the a quantum 
rotor model ||3^ that includes p interactions. The relation between the critical properties 
of the quantum versions of p spin-glass models and the experiments in pSf has been put 
forward in In addition, the connection between the static calculation supplemented 
by the marginality condition and the analysis of the out of equilibrium dynamics in contact 
with an environment developed in was also discussed in |Q . However, the reason why 
the transition changes from second to first order close to the quantum critical point was 
not clear from this analysis. It is one of the aims of this article to clarify this point, and 
study to what extent one can claim it to be general, with the use of the TAP approach. 
Quantum TAP equations for the SK model in a transverse field have been presented by 



Ishii and Yamamoto [R3] and Cesare et al [35|. The former use a perturbative expansion 



of the free-energy in the strength of the transverse field, and then follow closely TAP's 
techniques; the latter implement a cavity method. The TAP equations derived by Rehker 
and Oppermann [T^ for the fermionic spin-glass model coincide with the ones presented 
by Yokota for the classical GS model since these two models are thermodynamically 
equivalent ^ . 

Hence, our aim is twofold. On the one hand we want to present a quantum extension 
of the TAP approach to the statistical properties of disordered systems. Thus, after a 
short revision of the classical TAP approach in Section |^, we discuss in Section ^ the 
derivation of the quantum TAP free-energy and TAP equations using a general approach 
that extends the ones developed by Plefka P and Georges and Yedidia |P . The advantage 
with respect to previous derivations of quantum TAP equations is that this method can 
be applied to any quantum disordered model and it allows to obtain the TAP equations 
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as well as the TAP free energy. In Section ^ we present, as an example, the TAP free- 
energy and TAP-equations for the quantum extension of the p spin spherical spin-glass 
model studied in [^, We show that the TAP equations can be easily related to 
the equations for the order parameter in the Matsubara replica approach and also to 
some of the equations appearing in the real-time dynamic approach. The TAP analysis 
of this model furnishes a benchmark to study the generalization to the quantum case of 
the methods and interpretations developed for classical systems. Section ^ is devoted to 
the second aim of this article. Via the TAP approach we show that the same type of 
phase diagram naturally emerges for all systems having a discontinuous phase transition 
in their classical limit (these are models solved by a one-step replica symmetry breaking 
Ansatz within the replica analysis). In particular we relate the first-order transition close 
to the quantum critical point to the structure of metastable states. Finally, we present 
our conclusions in Section |6|. 



2 The classical TAP equations: a short revision 

In this section we present a short revision of the classical TAP approach to mean-field 
disordered spin models. The classical TAP free-energy is the Legendre transform of 
the free-energy with respect to local magnetic fields, 

- f3F{f3,mi) = TTexp(^-(3H -Y,hi{si-mi)^ , (1) 

where hi are Lagrange parameters enforcing the condition (sj) = mj. —j3F{j3,mi) is an 
effective potential that depends on the local magnetizations. The Lagrange conditions 
—df3F/dmi = hi, called TAP equations, fix the local magnetizations as functions of the 
local magnetic fields [|. The solutions {m"} of the TAP equations are stationary points 
of F{(3,mi). If they are also stable (all the corresponding eigenvalues of the free-energy 
Hessian are positive), they are identified ^ with pure states, also called TAP states. 
This interpretation was put forward by De Dominicis and Young ||^ who showed that the 
partition function in the classical SK model can be written as a weighted sum over the 
stable solutions of the TAP equations: 

Z = ^exp(-/3F(/3,m")), (2) 

a 

where the index a labels different TAP states, m° is an A^-vector encoding the local 
magnetization in the solution a, F is the extensive TAP free-energy of such solution and 
the sum runs over all TAP solutions. Consequently, the static average of any observable 
can be computed from Eq. (^. At low temperatures the TAP free-energy has a large 
number of minima. If one groups different TAP states with the same free-energy in sets 
C then the partition function can be written as 

Z = J2W,[3)exp{-f3Nf) (3) 
c 

^Note that within the TAP approach one does not average over disorder from the beginning as in the 
rephea method. Consequently the TAP free energy depends on the particular realization of disorder. 
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where the factor A/'(/3, /) is the number of solutions with TAP free-energy density 
F[f3,m°')/N = f. One can now replace the sum by an integral and exponentiate the 
factor A/'(/?, /); this yields 

lim i-lnZ= hm i- In / rf/exp(-iV(/?/ - a(A /))) (4) 

where we have taken the continuous limit and introduced the complexity 

a{(3,f)^hm lln(Ar(/?,/)) . (5) 

The configurations that dominate the sum are those having a free-energy density such that 
it minimizes (3f — cr{j3, /). The identity between the partition function and the weighted 
sum over TAP solutions has been demonstrated for many others models [1T2| , ^ and it 
is generally believed to hold for any mean-field disordered system. 



In the following we focus on "discontinuous glassy systems" [|2^ that are character- 
ized by having a discontinuous transition (the Edwards- Anderson order parameter, ^ea, 
jumps) that is still of second order thermodynamically. Within the replica analysis of the 
partition function these models are characterized by a one-step replica symmetry breaking 
solution below a static transition temperature Tg and a replica symmetric (RS) solution, 
that corresponds to the paramagnetic phase, at T > T,,. However, for intermediate tem- 
peratures Tg < T < Td there are an exponential in number of non-trivial TAP solutions 
that combine themselves in such a way that the sum (HI) is identical to the RS result. 



The relationship between metastable states and replicas has been put forward in |H3 



1^, E9i. Indeed, consider x different identical systems ("clones") coupled by an attractive. 



infinitesimal (but extensive) interaction. When there exist many pure states all the clones 
fall into the same state and the free energy for the system of x clones reads: 

jim ^InZ. = jim ^In / rf/exp(-Ar(/5x/ - a(A/))) • (6) 

On the other hand the computation of the left hand-side of (^ can be performed within 
the replica formalism: 

-1 -1 -1 

lim In = lim In Z^. = lim In Zj} (7) 

N-^oc (3N N^oo f3N N^oc,n~^o (3Nn ^ ^ ^ 

where the overline represents the average over disorder. Since the attractive coupling 
between the x clones is infinitesimal, the computation of the right-hand side of (^ reduces 
simply to the calculation of limn'^o{x/n') In Z"', where the replica symmetry between the 
n groups of x-replicas {n' = nx) is explicitly broken. When the system is in the replica 
symmetric phase {Tg < T), this reduces to study one-step solutions non-optimized with 
respect to x\ 

~ }}^ W^T^^ f dfew{-N{f3xf - a{(3, /))) = xExtTg^J^^piq^^; (3, x) (8) 

where /rep is the free-energy computed by using replicas, ^ea is the Edwards-Anderson 
parameter and x is the breakpoint (or the size of the blocks in the replica matrix). For 
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simplicity we consider that the inter-state overlap go equals zero. The definitions of these 
parameters are standard in the replica approach P] and they will appear in the analysis 
of the quantum p spin-glass model in Section 3. Since the integral on the left-hand-side of 
(H) is dominated by a saddle point contribution, one finds that, for a given temperature, 
fixing the value of x is equivalent to summing over states with a given energy density /. 
The relationship between / and x reads 

(9) 

Note that within this framework one does not have to optimize with respect to x. Instead, 
x is a free parameter and, by changing the value of x, one can consider different groups 
of metastable states. 

The analysis of the TAP equations reveals three temperature regimes for discontinuous 
glassy systems: 

• High temperatures < T. The system is in the paramagnetic phase, the paramag- 
netic TAP solution rrii = 0, for all i, dominates the sum and fpm = — In Z/ {f3N). 

is the dynamic critical temperature. Above Td the dynamics starting from a random 
initial condition converges asymptotically to the paramagnetic solution. 

• Intermediate temperatures Tg < T < T^. The replica analysis of the partition 
function indicates that the system is still in the paramagnetic phase. However, the 
study of the TAP equations and the dynamics show that at Td the paramagnetic 
solution is fractured in an exponentially large in number of minima of the TAP 



free-energy |TI], |^, Indeed one can recover these results also by the replica 

method. A careful replica analysis shows that there exist one-step solutions in these 
temperature regime other than the paramagnetic one. These solutions are in one to 
one correspondence with groups of states with a given free-energy density (through 
the relationship @)). For instance, one can follow the evolution of the threshold 
states (the states with highest free-energy) by tuning the parameter x. For these 
states, X = 1 when T = Td and decreases at lower temperatures. Moreover, the 
dominant contribution to Eq. (^ is given by the states characterized by x = 1, i.e. 
those with free-energy density such that 

P^'-^. (10) 

of 

These are the threshold states at Td and other groups when T < Td. Hence, between 
Tg < T < Td saddle-point solutions (corresponding to x = 1) that are not absolute 
minima of F dominate the integral since their number scales exponentially with A^. 
The final result for the free-energy density in this temperature range coincides with 
the one of the prolongation of the paramagnetic solution (that actually does not 
exist!). A naive replica computation fails to signal the difference between a true 
paramagnetic solution and the ensemble of non trivial TAP solutions with rrii ^ 0. 
The dynamic approach detects the change in free-energy landscape at Td since the 
system cannot reach equilibrium for any temperature below Td 
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• Low temperatures T < Tg. At the static transition temperature the complexity of 
the TAP solutions, which dominates the sum (0), vanishes. The static transition 
appears as an entropy crisis since the part of the total entropy that is related to the 
large number of states disappears. For T < Tg the TAP states which dominate the 
integral in Eq. (H) correspond to the equilibrium glassy phase. Dynamically, does 
not play any role. The out of equilibrium dynamics is dominated by the threshold 
states, which are the highest ones in free-energy and that are characterized by flat 
directions in the free-energy landscape. 

Note that via the TAP approach one can obtain a reasonable justification of the 



marginality condition [501 often used to obtain information about the out of equilibrium 



dynamics starting from a pure equilibrium computation [^]. Indeed the value of x fixed 
by the marginality condition corresponds to the TAP states which are marginally stable 
(the threshold states): the flatness of the free-energy landscape around these states is 
responsible for aging p2| . 



3 The quantum TAP equations 

In this section we present a simple procedure to derive TAP equations for generic com- 
pletely connected quantum systems. We also expose the physical meaning of quantum 
TAP equations by the cavity method [Q. 

We are aware of two publications where TAP equations for quantum systems have been 
already presented ^ . With respect to these works our derivation is more systematic, 
simple and it allows one to obtain the TAP equations as well as the TAP free-energy for 
any completely connected quantum disordered system. 



3.1 Formalism, notations and models 

The formalism that we use to derive TAP equations for generic quantum problems is very 



similar to the one described in g y] and, it follows even more closely, the one used in p4 
to obtain the dynamical TAP equations for classical disordered models. 
We focus on systems characterized by the potential energy: 

Hp = - E J2 Jiu-,i.€ ■ ■ ■ 'I z = l,...,Ar a = l,...,m (11) 

ii<...<ip a 

where s, may represent an SU(2) spin (m = 3), a rotor (m > 1), a spherical spin (m = 1) 
or a space-coordinate (m = 1) and the couplings between the different Sj, are 

independent random variables with zero mean and variance 

As a consequence the following derivation applies to (completely connected) Heisenberg 
models, quantum rotor models and quantum continuous systems^. Without loss of gen- 

^ The results obtained in this section can be straightforwardly generalized to more complicate poten- 
tials containing different monomials or characterized by a more complicate tensorial coupling between 

e"l ... 
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erality and to simplify the notation we shall suppress the index a in the rest of this 
section. 

For classical spin-glasses TAP showed that all the local magnetizations rrii, i = 
1, . . . . ,N, are needed to derive the relevant free-energy density to describe the metastable 
properties If one is interested in the dynamics of classical disordered mean- field sys- 
tems, one has to Legendre transform not only with respect to all time-dependent local 
magnetizations mj(t), but also with respect to the autocorrelation C(t, t') = l/NJ2i{si{t)si{t')) 
and the linear response R{t,t') = l/NY.^5{si{t))/5h,{t')\h=o 0. 

In order to describe the metastable properties of a quantum disordered model we 
shall show that it is necessary to Legendre transform with respect to the local average 
coordinates, mj(r), and the autocorrelation function in imaginary time, C(r, r'). The 
quantum TAP free-energy reads 

-/3F(/5,m,(r),C(r,r'),«)L=i = 

\nJVs{T)exp --J dr {Hk + aHp{s)) 



2h- 

+ 



-2 / dr dr' EA(r,r') (C(r, r') - ..(r)s,(r')) 
n Jo Jo ~r 



1 [^^ 

rfr}J/ii(r)(mi(r) - Si(r)) 



(13) 



where is the kinetic energy, X>s(r) indicates the functional measure on the configuration 
space and a is a parameter whose role will be clarified in the following. For instance, if s, 
are SU(2) spins if^ is the Berry phase and the functional measure is restricted to periodic 
functions Sj(r) (with period /5) satisfying the constraint s^(r) = 1. The sources hi{j) and 
A(r, r') have the role of Lagrange multipliers fixing the average value of the coordinates 
and the correlation: 

m.(r) = (s.(r)), (14) 

i 

Once the TAP free-energy F is known, one can derive the TAP equations as Legendre 
relations, 

hf\ 2 6(3F 



5mi{T) " N5C{t,t') 

Until now we have not not used the scaling (|T2]) and all these definitions can be equally 
applied to finite dimensional systems. The great simplification due to the mean field 
character of the interactions in ([l2[) is unveiled if one performs a perturbative expansion 
of Eq. (|1^) in a and writes —PF{P, mi^r), C{t, t'), a) as a power series in a: 



n^.n i ,r( '\ . ^1 9"(-/5F(/5,m,(r),C(r,r'),«)) 
/5F(/3, mi(r), C(r, r ), a) = 2^ 



n=o ^' '^a" 



a" . (17) 

0=0 



In fully-connected models, if one chooses the correct order parameters (which are mj(r), 
C{t,t') in the quantum case), the perturbative expansion (^) around the pure kinetic 
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theory is actually a simple sum over three terms. Higher order terms in the series vanish in 
the thermodynamic limit due to the scaling of Jji,...,jj, with respect to A^. In more general 
cases, in finite dimensions, this will not be the case and ([T^) becomes a l/c? expansion 
around mean field theory, where d is the spatial dimension |Q. 

Let us consider in more detail the terms arising from the expansion (plTl). The zeroth- 
order one is simply — mj(r), C(r, r'), 0), i.e. the free-energy of free spins con- 
strained to have local magnetizations mi(r) and a global correlation function C(r, r'). 
This term depends only on the nature of the degrees of freedom, whether they are SU (2) 
spins, rotors or space-coordinates. In particular it can be analytically computed only 
if Si are spherical spins or space-coordinates. In the other cases one has to resort to 
approximations or numerical computations. 

The first-order term is the naive mean field free-energy: 



d{-(3F) 



da 



a=Q ^ ii<...<ip 



■ ■ mi 



i\<...<ip 



Note that the decoupling of the spins for a = is essential to obtain the last identity. The 
second-order term depends on the correlation function and the overlap function Q(r, r') = 
'mi{T)mi{T') / N only and equals 



dr dr' (C^(r, r') - g^(r, r') - p(C(r, r') - g(r, t'))Q^'~\t, r')) , (19) 



Ah^ Jo 

Using the scaling of the couplings with and by the same arguments developed for 
classical systems 0, we have verified that all orders n > 3 in the series (|T7p are 
suppressed in the thermodynamic limit. 



3.2 A cavity interpretation 

First of all, let us write the TAP equations in a way which allows one to clarify the 
physical meaning of the different terms: 



S{-PF) 



5mi{T) 



2 Si-PF) 



a=Q 



1 fP^ 



h Jo 



dr' 



p{p — 1 



{QiT,T')-C{T,T'))Q^-\T,T') 



(20) 



N 6C{t, t') 



Qcav^^^ r') = I [gP-^(r, r') - C^^'ir, r') 



a=0 



(21) 



The solutions to these equations are expected to be time-translation invariant since we 
are developing a description of equilibrium and metastable properties. Therefore hf"''" 
is indeed independent of the imaginary time and G'^"''" depends only on the difference 
between r and r'. 

An understanding of the meaning of the quantum TAP equations follows from the 
analysis of F for a = 0. Indeed, by tracing out all the spins except Si in the partition 
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N-1 spins + cavity 



cavity + Weiss fields 




Figure 1: A schematic representation of the action of the — 1 spins on the cavity which, 
due to the infinite connectivity, simply reduces to hi"""" and G™''. 

function produces a single-site measure (for si) whose action reads: 

- \ [H,{s,{t)) + hr{T)s,{T)] - ^ ^ dr' s,{t)G^%t,t')s,{t') , (22) 

n Jo 2h Jo Jo 

where Hki^Siij)) is the kinetic energy for the spin Sj. As a consequence the TAP solutions 
are the self-consistent relations that relate G'^"'^{t,t') and /i™"(r) (which are functions of 
C{t,t') and mj(r)) to C(r, r') and miij) obtained from the single-site action (p^). 

Equations ( pO|) and (|2l|) show that the action on the zth spin of the A^ — 1 remaining 
ones reduces simply to hi"""" and G™''. This implies that tracing out all the spins but 
the zth one produces a Gaussian measure for the instantaneous magnetic fields hiij) = 
~J2i2<...<ip-Jii,..;ip^h'''^ipy whosc mean and connected two-point correlation function 
equal respectively h'^'^^ and G^°''"{t — t'). 

The expression of hf^^ and G'^°''"{t — r') can be justified within the cavity method [Q. 
Let us focus for simplicity on the p = 2 case for which 

hr = -T.J^Ask)N-i (23) 

k 

where {■)n-i represents the thermal average with respect to the system with the zth site 
removed. {sk)N~i is not simply equal to ruk, which is the mean magnetization for the 
system of N spins. A correction term, first discovered by Onsager, appears: 

{sk)N~-i = mk- T / dT———Ji,kmi{T') = mfc - - / dr [C(r) - Q] Ji^km (24) 

II Jo l>,k\J ) ' "'O 

Plugging Eq. (p^ into Eq. (p3|) and using the scaling of the couplings with A^ one recovers 
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the expression for /i^"" given in (^0|) in the p = 2 case. Whereas for C^"""" a similar 
computation ^ gives back the expression given in (pi]). 

Finally, we remark that the main difference between the classical and the quantum 
TAP approach is that in the latter the cavity interaction consists not only in a cavity field 
but also in the "Weiss function" G'^°'^{t — r'), which is a function of (imaginary) time. 
This already happens in the mean-field theory of quantum non-disordered systems p2[ for 
which local quantum fluctuations are taken into account exactly, whereas the spatial ones 
are frozen. For disordered systems, even in the limit of infinite dimensions, one has to take 
into account not only the local quantum fluctuations but also some spatial fluctuations: 
all the instantaneous magnetic fields have the same variance but their averaged values 
fluctuate from site to site. 



4 A continuous disordered quantum model 

In this Section we apply the method of Section ^ to the study of the quantum spherical 
p spin-glass model. We derive and analyze the TAP free-energy density and the TAP 
equations for the local magnetization and correlation function in imaginary time. We 
relate these equations to the equation for the order parameter in the Matsubara replicated 
approach to equilibrium and in the Schwinger-Keldysh approach to the non-equilibrium 
dynamics. 



4.1 The model and its TAP equations 

A model of a quantum particle with position s and momentum p that moves on an 
iV-dimensional random environment is defined as 

p2 TV 

H[P,S,J] = —- Jh-i,Sn-Si, ■ (25) 

h<...<ip 



A Lagrange multiplier z enforces the averaged spherical constraint 

^ E(^?) = 1 • (26) 



1=1 



The random interaction strengths Ji^,„ij, are taken with zero mean and variance defined 
in Eq. (|r^). This model is a possible quantum extension of the spherical p spin-glass 



model introduced in |^ and it is a particular realization of the class defined in (|TJ 
corresponding to space-coordinates Si constrained to move on a N-dimensional sphere. 

The zero-th order term of the expansion (|17D can be readily computed for this model. 
By setting a = 1, rescaling time according to t ^ rh/ J, and defining the "quantum 
parameter" F = fi^/(JM) we obtain the following expression for the quantum TAP free- 
energy ([T|) : 

-/5F = -Trln(C-Q) + -Tr(^^j + dr ^n.-v'^nM • • • m,^(r) 



^Indeed easier since there is no reaction term for G"^"" on a completely connected lattice. 
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+ J I dr dr' (C^(r, r') - Q^{r, r') - p{C{r, r') - Q{r, r'))Q'-\T, r') 

dTz{T){C{T,T) - 1) 



'0 JO 

NI3 [1^ 



2 Jo 



(27) 



The physical parameters mi{T) and C(r, r') are fixed by the quantum TAP equations (|TBp 

hi{T)= Ji,i2,...,ipT^i2i'r) ■ ■ ■^ipi'r) (28) 

i2<...<ip 



+ / dr' 



-{C - Q)-\t, t') + ^^^-^ (Q(r, r') - C(r, r')) Q'-\r, r') 



z{r)6{r - r') = {C - Q)-\r, r') + 5{r - r')^|^ + | 



mi(r') 



Finally, setting /ij(r) = and using that at stationarity, mi{T) = rrii, Q{t,t') = ^ea, 
z{t) = z and C(r, r') = C(r — r'), the previous equations are simplified to 



1 d'^Cir) 
f 9r2 



"2 Jo 



/3 



rf^' (C^-i(r-r')-gg-^) (C7(r')-gKA) 
+;.(C(r)-gEA)-5(r), 



(29) 



ZTfl^ ^ ^ Ji,i2,...,ip'^i2 ■ ■ ■ '^ip 

j2<---<jp 

rn. I rfr (C^-i(r') + (p - 2)^^"^ - (p - l)C(r')gL-') • (30) 



4.2 Analysis of the quantum TAP equations 

In the classical case the TAP equations admit a large number of solutions at low tem- 
peratures. In the following we shall show that this remains the case in a certain regime 
of T and F. Furthermore we shall classify them by their Edwards- Anderson parameters. 
Finally, we shall exhibit several properties of the TAP solutions valid at low temperatures. 



4.2.1 A simple equation on the Edwards- Anderson parameter 

Let us analyze in detail the equations for the local magnetizations, mj. First of all we 
note that a simple equation that relates ^ea to the potential energy density derives from 
Eq. (p8|). In fact, by multiplying Eq. ( ^8]) by mi/N and summing over i = 1, . . . N one 
obtains (for hi = 0) 

= - ^,,,^'^ +^ E ^.....,m,...m,^-^^^(C'(0)-/5gEA)gL-' (31) 
C(0) - /3gEA ^ n<...<v 2 V ^ 

where we introduced the discrete Fourier transform of the correlation 

C(uj) = f dre'^^C^r) . (32) 
Jo 
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Following Kurchan et al we introduce the an^ftt/ar variables = rrii/ ^Jq^ and define 
the angular potential energy density 

= Jh-i^<^h---(Jip ■ (33) 

n<...<«p 

For any fixed energy level S, Eq. (|3TD becomes a second-order polynomial equation for 
^ea; the solution is determined by 



q'J'-\CiO) - Pq,,) =z^ = --^{ -£{a) ± ^ S\a) - ) (34) 
and £^TH, called the threshold value, is given by 



= -f-^ . (35) 

The right-hand-side of Eq. ( p^ has to be real. This imposes the condition £ < £^th since 
is a negative quantity. 
For each sign in Eq. (pi]), its left-hand-side has a bell-shape as a function of ^ea- It 
vanishes at ^ea = and ^ea = C{Qi)/ {3 and attains its maximum at ^ea = (1 — 2/p)C(0)//?. 
Hence, at fixed values of £ and T, Eq. ( |34D has none or two solutions, ^ea = with 

0< g' <^ ^^^^ 

(1 - 2/^)^(0) C'(O) 

/3 - ^ - /3 ^^^^ 

(we assume, as expected, that Cir) is positive for all r). In the classical case, the minus 
sign in Eq. (|3^) leads to a value of ^ea that is a minimum of the TAP free-energy for 
all £ < £^TH- The Edwards- Anderson parameter determined in this way has the expected 
physical behavior [12|. In Appendix A we show that in the quantum case one has to 
choose the minus sign in Eq. (0), too. Thus, ^ea is determined by 



Q'l''\CiO) - /?gEA) = ^ y-E{o) - ^£\a) - £^ . (38) 

This equation still has two solutions. It can be proven that the solution with the larger 
absolute value of ^ea has the correct physical properties. In particular, it is connected to 
the classical solution, and it is then the solution to be kept. Thus there is a one to one 
correspondence between ^ea and £. 

It is of particular interest, as we shall show below, the threshold solution £ = £^th- In 
this case the equation for ^ea becomes 

1 = (C'(O) - /3gEA)^gL-' • (39) 

Note that this equation coincides with the one found with the Matsubara formalism 
using the marginality condition to fix the block size x in the replica matrix |38|. Fur- 
thermore, it coincides with the equation for the dynamic value of the Edwards- Anderson 
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parameter ^ea = linit^oo linit„^oo + tw,tw) obtained from the study of the real-time 
dynamics of the quantum model evolving in contact with an Ohmic quantum environment 
p7| , when one takes first the thermodynamic limit, next the long-time limit of the sys- 



tem's dynamics in contact with the environment and, finally, the strength of the coupling 
to the environment to zero. The relationship between TAP, Matsubara and dynamical 
approach will be discussed in Section I 



4.2.2 Multiplicity of TAP solutions 



The equations ( p4D reveal an interesting structure of the TAP equations ( poD and (pO|). 
For a given value of the angular potential energy £ the TAP equations decouple in two 



different sets: Eqs. (^) and (|3^), with the spherical condition on C, determine the corre- 
lation function, the spherical parameter and the Edwards-Anderson parameter; whereas 
Eqs. (|30D determine the angular variables only. They read 

= -pS{(T)ai=p Ji,i2,-,ip(^i2---(^ip ^ (40) 

i2<---<ip 

, - ---^^^<.- + ^<5(o),'--E(o). (41) 

where we have defined 



P 



S(0) = |y^ C^-\r). (42) 



For a given value of the angular potential energy £ Eqs. (^0]) allow one to determine the 
angular part of the TAP solutions. 

In general for a given value of £^ Eqs. (^) and (^) determine the correlation function, 
the spherical parameter and the Edwards-Anderson parameter in a unique way. (An 
exception to this rule are the paramagnetic solutions which however do not correspond 
to any £.) As a consequence, the multiplicity of TAP solutions is entirely due to Eq. ( ^OD 
which, for certain values of can admit an exponential (in A^) number of solutions Af{£). 
The complexity function of £ is then defined as 

1 



ai£) ^ hm -MAfi£)) , (43) 
Equation (^) already appears at the classical level. The complexity has been computed 



by Crisanti and Sommers [jl2[ and Cavagna et al. [|T^ with the following result. There 



are typically no solutions for £ < £eq, whereas for £eq < £ < £th the complexity reads 



TH 



for £eq<£ < £t:h (44) 

where £eq is the value at which (j{£) vanishes. A plot of this function is traced in Fig. ^ 
for p = 3. 
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Figure 2: The complexity (j[S) as a function of S in the interval S^q < S < S-yu for p = 3. 



4.2.3 A low temperature and low F approximation 

In the classical case the TAP equations separate in two sets: N equations for the angular 
variables and one for the Edwards- Anderson parameter. The former admit an exponential 
number of solution and are studied from a statistical point of view (one computes the 
number of solution, and the typical properties of solution corresponding to a given S), 
whereas the latter can be easily solved. In the quantum case the analysis of the equations 
for the angular variables is identical to the one used for classical systems. The analog of 
the equation that determines ^ea becomes now a differential equation for C (r) that has to 
be studied numerically. This differential equation can be mapped exactly onto the ones 
analysed with the replica method, as we shall show in Section ( [4 .31 ), and its numerical 
solution can be found in ||3^. In this section we perform a low-temperature and low F 
approximation, also discussed in p8[, that allows one to obtain some qualitative results 



that remain valid for the exact solution. 

At low temperature and low F, the extension of the imaginary time-interval diverges 
[0,/3 oo] and the periodic correlation C(r) is expected to have a rapid decay, over 
a short time-interval, from 1 to its "asymptotic" value, say, at r = (3/2. Moreover the 
"regular" part of the correlation, that we define as 



(45) 



can be assumed to be small. Therefore we can expand the TAP free-energy in powers of 
Qreg{t)- Up to terms of the order of qREciTY we obtain 



(3F 1 



N 



IYln(gj,EG(r)) 



2F I 9r2 



\- T. 

ii<...<ip 



rrii 
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J' drql,,{T) - L (q^UO) + Q.. - 1) , (46) 

where we have focused on the space of time translation invariant (TTI) functions (since 
the TAP solutions are TTI this does not imply a loss of generality). 

Within this approximation the TAP equations become quadratic in Fourier space, 



+ q..oM + qlf qLcM = , (47) 



and yield 



z + ul/V ± d{z + ullV'f - 2p{p - l)gi 



,p-2 
EA 



q^EcX^k) = — -, —jr2 • (4^ 



By taking Uk = and comparing to Eq. (|38D one obtains 



l-p/2 

S = -^S^ . (49) 
P 

The spherical constraint reads 

1 - ?EA = \ ^qREG{^k) = [ — x'iuj) coth ( ^ ) (50) 



where 



l-p/2 

X"{uj) = Im q^Eci^k = -i^) 



2 l-p/2 

+ . (51) 



p-i\ \ 

The integral in Eq. (|50D has to be taken on the interval ui G cu+J such that the square 
root is real. 

In the low temperature limit, we approximate coth(/3c<j/2) ~ 1 and, by changing 
variables in the integral, we obtain 

m£,p) = !^!i^^(l - qM.-'^^' (52) 

p 

with 



I{S, p)=2 dx V^t'h - (^ + ■ (53) 

Equation (0) yields a relation between F, ^ea and S of the form 

ri\£,p) = ct{l-q,,fq^Z-'^/\ (54) 

with ct a numerical constant. For each S, there is a solution with a physically meaningful 
value of ^ea that is close to 1, until reaching a critical Fmax(^^)- This value tells us when 
the TAP solutions associated to £ disappear. It can be easily proven that I{S,p) is a 
growing function of S; hence, Fmax(^^) is a decreasing function of S. This implies that the 
TAP solutions that are at the threshold level disappear first than those that are at lower 
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Figure 3: TAP free-energy versus F at zero temperature. The curve (1) corresponds to 
threshold states which are the first ones to disappear at F^h; the curve (3) corresponds 
to the states with the lowest free-energy which are the last ones to disappear at Frsb and 
the curve (2) corresponds to an intermediate state. 



values of £. This is again similar to the dependence of the classical TAP solutions with 



temperature |Tl|: the solutions corresponding to the threshold level disappear at a lower 
temperature than the ones corresponding to the equilibrium level Tmax(^^th) < ^max(^eq) 
and, more generally, Tmax(^^i) < Tmax(^2) if £i > ^2- 

The low frequency behavior of the spectral density x"{^) = QnEai^k = —iuj) of the 
threshold states is gapless, 

X'iuj) r^to for w ^ 0+ , (55) 
whereas all the other states [S < S-yh) have a gap A in their excitation spectrum. 



x'iuj) ~ Vuj-A for ^ A+ . (56) 

Furthermore we have studied the dependence of the free-energy on S and F in the 
low temperature limit. Plugging the solution (^8|) into (|46|) we find, after a tedious 
computation. 



N 




2sinh — cux"(cu) + \ p-2-^\8 (57) 




where ^ea satisfies Eq. ([38|). This expression allows one to study the evolution of the 
free-energy of the TAP states as a function of F. We have found that if one knows a TAP 
solution at zero temperature and zero F one can follow it continuously in F. As in the 
classical case TAP solutions do not cross, merge nor divide in this model. 

Figure ^ summarizes these results in a schematic way. Finally, note the special role of 
the threshold states, which are gapless (contrary to the others), the first ones to disappears 
and the ones with highest free-energy density. 
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4.2.4 Stability of TAP states 

In the classical case one can check the stability of TAP states. In the quantum case this 
is a difficult task that has to be performed numerically. In the following we shall limit 
ourself to prove that the TAP solutions characterized by £^ = £^th (threshold states) are 
characterized by zero modes and are hence marginally stable. We expect that a complete 
stability analysis will con&m that the TAP states characterized by £^ < S^u are stable. 

Let us focus on the reduced free-energy Hessian d"^ F / dmi{T)dmj{T') evaluated in a 
TAP solution {mf }. This matrix depends on r, r' only through their difference. Therefore 
it is diagonal in Fourier space. Focusing on zero frequency, the original problem reduces 
to the diagonalization of the following matrix: 

^i,3 = - Jw3,-,ipKs ■ ■ ■ Kp - p^(iEA~^kj > (58) 

i3<...<ip 



where ^ea = HA^tY I ^ ■ The density of eigenvalues of A has been computed in and 



except for the isolated eigenvalue corresponding to the eigenvector mf , it is a semicircular 
law centered in —p£(^^^ ^ with width —pE^^q^^^ ^ . Consequently, threshold states are 
characterized by a vanishing fraction of zero modes. 

4.2.5 The classical limit 

The classical limit of Eqs. (|29|) and ( pOD yields the classical TAP equations computed by 
Kurchan et al. In fact, in the classical limit, C(r) = 1 and the parameter z is fixed 
by integrating Eq. (|29D between O"*" and (3^ . This yields 

p(3 



(3{l - q) 2 
By inserting this value of z in Eq. ( pO]) we obtain 



+ f (l - f-') ■ (59) 



+ - q)q^ 2^ mi = p ^ Ji,i2,...,ipmi^ ■ ■ .rriip , (60) 

that coincide with the classical TAP equations for the local magnetizations. The equation 
that fixes ^ea as a function of £ in the classical limit is simply obtained from Eq. (0) by 
setting (7(0) = 1. 

4.3 Relation between TAP, Matsubara and dynamic approaches 

In Section |^ we have recalled the relationship between TAP, replica and dynamical ap- 
proaches in the classical case. In this Subsection we show how these connections are 
generalized to quantum systems. 

4.3.1 TAP-Matsubara 

Via the replica analysis in the Matsubara imaginary-time framework and within a Istep 
RSB Ansatz, the order parameter is the n x n matrix Qab which is fully described by: 



19 



n identical diagonal elements ^'^(t) that depend on the the imaginary time r, n{x^ — 1) 
constant elements ^ea that occupy the x x x blocks around the diagonal, the remaining 
71^ — n{x'^ — 1) elements go that in the absence of an external field are identically zero. 
In the — > limit, the three parameters gd(r), ^ea and x, together with the value of 
the Lagrange multiplier that enforces the averaged spherical constraint determine the full 



solution of the problem |38 



The connection between TAP and the Matsubara approaches is obtained by identifying 
the Edwards- Anderson parameters ^ea in the two approaches, C (r) with the r-dependent 
diagonal parameter qdij) in Qab, and the Lagrange multipliers. In particular we have 
shown that subtracting the equation obtained for a ^ h from the one corresponding to 
a = h one obtains Eq. (|29|). In the Matsubara approach one has another equation for ^ea 
in which x acts as an external parameter. Therefore by fixing the value of the breakpoint 
one fixes the value of ^ea- As in the classical case two different recipes to fix the breaking 
point parameter x, namely optimization and the marginality condition, lead to the static 
and dynamic transitions, respectively. In the TAP approach the role of x is played by 
£ that enters the equation for gEA as a parameter. We have found that the relationship 
between x and £ is encoded in 

^ ^ . (61) 

of 

where a is the complexity defined in (H), see Appendix B. This suggests that in the a 
quantum problem the relationship @ is generalized to 

_ lim ^Ve-^^^^" = - hm -^In / d/e^(-^-^+'^(^'^)) (62) 
N^oo (3N ^ N^oo (3N J •' ^ ^ 

= 2;Extrgi^^ g^(T-)/REp(gEA) 

Using Eq. (|6TD we have found that the Matsubara equations for ^ea and ^^(t) and the 
TAP equations for ^ea and C(r) coincide. For instance, the TAP equations for the highest 
TAP states (threshold states) and the lowest TAP states coincide with the ones obtained 
in the Matsubara approach by using the marginality condition and the extremization 
with respect to x, respectively. Moreover, as another confirmation of Eq. (|62D we have 
verified that the free-energy obtained from the Matsubara computation equals the 
one obtained in the TAP approach for all values of (3 and P. In other words, we have 
checked that 

- /3F = In Y.^''^^" ■ (63) 

a 

As a consequence the phase diagram that follows from the TAP approach coincides with 
the one obtained in 



4.3.2 TAP-out of equilibrium dynamics 

The study of the real-time dynamics of the p-spin quantum model evolving in contact 



with an Ohmic quantum environment has been performed in |^ . The dynamical behavior 
is characterized by two regimes. At high temperature and high P the system equilibrates 
in the paramagnetic state via an equilibrium djTiamics. Whereas at low temperature and 
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low r the systems ages and remains out of equilibrium also at infinite times. As in the 
classical case we have found that the long-time out of equilibrium dynamics is dominated 
by the threshold states. This is proven by the fact that the equations for and C(r) 
within the TAP approach coincide with the ones obtained from the study of the real-time 
dynamics, when the following limits are taken in its precise order: lim^^o liint^oo liniAr^oo- 
In order words, when one takes first the thermodynamic limit, next the long-time limit 
of the system's dynamics in contact with the environment and, finally, the strength of 
the coupling to the environment 7 to zero. Notice that this equivalence holds for the 
paramagnetic states also. Finally, we have shown that the relationship between the ef- 
fective temperature arising in the asymptotic out of equilibrium regime and the 
complexity is, as classically, 

1 da{(3,f) 



df 



(64) 



/=/t, 



Note that the connection between TAP and real-time dynamics is done by identifica- 
tion of several equations. A more precise analysis, along the lines of [Q, should prove 
the full equivalence of the two methods. 



5 The phase diagram of discontinuous glassy systems 

In this section we present some general arguments that allow one to predict the phase 
diagram of discontinuous glassy systems. Since the free-energy landscape plays a key role, 
we expect these results to have a certain degree of universality and to apply to this entire 
class of disordered systems. 

5.1 The static transition 

Let us focus on two limiting regions of the phase diagram: around the classical phase 
transition (F = 0) and around the quantum phase transition (T = 0). 

In the former case the physics is well known and it is reviewed in Section |[ The 
effect of switching on weak quantum fluctuations consists only in a weak variation of 
the complexity a. For this reason the effect of quantum fluctuations reduces simply to 
a variation of the thermodynamic (T^) and the dynamic (T^) transition temperatures 
(respectively lines (1) and (3) in Fig. ||). 

At zero temperature and low F the system is in the glassy phase (GP), whereas at 
very high F quantum fluctuations destroy the glassy phase and the system is a quantum 
paramagnet (QPM). As a consequence one expects that a quantum phase transition should 
divide these two regimes at a certain value Fc. 

At zero temperature the complexity is expected to remain a smooth function of the 
free-energy density^. Consequently, equation (^) implies that the sum over the exponential 
number of glassy states is always dominated by the lowest ones in free-energy since the f5f 
term in the exponential largely dominates in this limit (this is different from the classical 
problem in which other states dominate between Tg and T^). At zero temperature these 

^For instance the complexity does not blow up for T — )■ since it is bounded by the logarithm of the 
number of energy minima divided by N, which is a finite quantity independent of temperature. 
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Figure 4: A schematic representation of the phase diagram which is expected to be generic 
for systems having a discontinuous transition in the classical limit. The line (1), from 
(Tc, 0) to the point B, represents the static transition between the classical paramagnet 
(CPM) and the classy phase (GP). The region in red between the lines (1) and (3) is the 
phase in which the CPM is fractured into an exponential number of TAP states. The 
line (2) from the point B to (PcO), signals the static transiton between the quantum 
paramagnet (QPM) and the GP. The line (3) indicates the dynamic transition as a 
function of P. 

are the states with lower angular potential energy, i.e. with S = E^q. Now, from Fig. ^ 
we conclude that o"(£^eq) = at zero temperature and for all P. For this reason the 
mechanism behind the transition must be totally different from the classical one. The 
transition cannot be related a configurational entropy that vanishes when approaching P^ 
from above ( "entropy crisis" ) since this quantity is always zero at zero temperature. 

Indeed, according to Eq. (|61|), if we assume that da{P, f)/df < +oo when T — > 0, 
then X — > for all P in the glassy phase. In the paramagnetic phase instead x = 1. Thus, 
X must jump at the transition. If the Edwards- Anderson parameter also jumps at Pc, the 
susceptibility is discontinuous, and the transition is of first order thermodynamically. As 
in the previous case the effect of switching on thermal fluctuations reduces simply, for low 
T, to a variation of Pc (line (2) in Fig. 

Another hint on the difference between classical and quantum phase transition can be 
gained by a technical remark. It is well known that the paramagnetic solution of the clas- 
sical problem remains stable in the low temperature phase. This is a spurious solution of 
the mean field equations which has to be discarded in the analysis of the low temperature 
regime. In the quantum case, one also expects to find a spurious paramagnetic solution, 
which is the continuation of the classical paramagnet to low temperatures. This solution 
exists to the left of line (1) in Fig. H, consequently one expects coexistence of two para- 
magnetic solutions: a physical one which is the continuation of the quantum paramagnet 
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valid at low temperatures and high T and a spurious one which is the continuation of the 
classical paramagnet. 

In the classical case the transition is of second order even if the order parameter jumps 
discontinuously. This peculiar behavior is due to the fact that near the transition the 
paramagnetic state is fractured into an exponential number of states which continuously 
become the ones responsible for the glassy phase at low temperature. This is not possible 
at zero temperature (the quantum paramagnet is not formed by a collection of glassy 
states) and therefore it is reasonable to expect a quantum first order phase transition 
between the glass phase and the quantum paramagnet .0 

Finally, note that Fqp^ = F^p = Fcpm on the point B. We then expect that a first order 
transition line separating the QPM from the CPM starts at this point. This line should 
end on a point C given that for very large value of F and T the quantum and thermal 
fluctuations are so strong that the system becomes non interacting and in this case only 
one paramagnetic phase exists. In the analysis of the p = 3 spherical spin-glass model 
[p8| the line BC has not been found. We conjecture that in this case the line BC is so 
short that it is very difficult to find numerically. Furthermore, within the accuracy of the 
algorithm, the dynamic and static critical lines collapse at the point B. In the quantum 
model studied in instead this line has been detected and it was demonstrated in this 
paper that its length increases with p. 



5.2 The dynamic transition 

Now that the equilibrium phase diagram is completely predicted from a qualitatively 
point of view, we can focus on the non-equilibrium regime. As noted previously, low 
quantum fiuctuations simply change the values of but do not change qualitatively the 
dynamic transition which remains second order in the sense that the asymptotic energy is 
continuous across the transition, but its derivative is not. This remains true until the line 
(2) reaches the line BC. After this point the dynamic transition between the the quantum 
paramagnet and the threshold states becomes first-order, i.e. the asymptotic energy is 
not continuous across the transition. This, of course, is very difficult to see numerically 
since the discontinuity has a very small value. 



5.3 Summary 



In summary, through some general arguments based on the TAP approach we have pre- 
dicted a phase diagram that should have a certain degree of universality since its form 
is determined by the qualitative form of the free-energy landscape. Indeed, not only 
the quantum p spin spherical model exhibits the phase diagram displayed in Fig. ^ but 
some other classical and quantum models share exactly the organization of phases and 
transitions 
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36, 45 



^Note however that other scenarios are possible. For example the number of glassy states could 
diminish when F increases and vanish exactly at Fc. In this case the glassy states could be grown up 
continuously from the quantum paramagnet. 
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6 Conclusion 



In this paper we have derived TAP equations for a large class of mean field disordered 
quantum system. Moreover we have applied the TAP approach to the quantum version 
of the spherical p spin model. The study of this system, whose real-time dynamics and 



statics have been analyzed in has furnished an ideal benchmark to generalize 

to the quantum case several concepts developed for classical disordered systems. Armed 
with this knowledge, founded on the study of the free-energy landscape, we have shown 
that the same phase diagram, presented in Fig. |^, naturally emerges in a large class of 
quantum disordered systems, those having a classical discontinuous transition. 

Whether other models like the SK model in a transverse field, or its soft spin version 
studied in also have such crossover in the transition from the disordered to the 



ordered phase is an issue that deserves revision. For the moment, no study of models 
with classical continuous transition has shown this feature. However, it might have been 
masked by the methods used in previous studies. The soft SK model might be the easiest 
example where to answer this equation via, e.g., a careful application of the rephcated 
Matsubara approach. 

We would like to stress that the TAP approach furnishes an alternative and more 
transparent route to replicas which has also the advantage of showing explicitly the weak- 
ness of the mean field description. Let us cite one example. The marginality prescription 
in the replica approach becomes the more transparent statement that the non-equilibrium 
dynamics is dominated by the TAP states which are marginally stable, i.e. the flatness of 
the free-energy landscape around these states is responsible for aging. Concerning one of 
the weaknesses of the mean field description we would like to underline that the enormous 
number of pure states (with different free-energy densities) found for mean field models 
cannot persist in finite dimensions and the majority of them should become metastable 
states. How this changes the mean field scenario is an active domain of research for 
classical systems |l54 |. 

We remark that interesting continuations of our work concern, on the one hand, the 



application of the TAP approach to different quantum mean-field models pi], E^, 



and, on the other hand, the generalization of the static quantum TAP approach to real- 
time dynamics (for classical systems this has been done in [24||). This would allow one to 



show the relationship between long-time dynamics and free-energy landscape for quantum 
systems directly. Finally, the precise definition of a "quantum state" is a delicate matter 
and merits further analysis. In this paper we have simply called "state" a minimum of 
the TAP free-energy density. One possible way to verify the existence and stability of 
these states is by studying the dynamics of this system starting from particular initial 
conditions as done in |HS|, for the classical model. This study is underway 
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Appendix A 



In tliis appendix we sliow tliat tlie equation for tliat leads to a solution with the 
correct physical properties is Eq. (|34D with the minus sign. Indeed, we search a solution 
that corresponds to a minimum of the TAP free-energy. The full stability analysis, that 
involves the evaluation of the complete Hessian of the TAP free-energy, is rather hard 
and has to be done numerically (for instance, the form of C{uj) can only be obtained 
numerically) . 

However, we can still perform a partial analysis that suffices to justify the choice of the 
minus sign. Let us concentrate on the following diagonal elements of the Hessian: 



di-PF) 



1 



pip 



(C(0) - /3gKA)2gEA 



•iEA 



-1 



dC{0)dC{0) 
From 2;-|-'s definition we find 



{C{0)-qf3y 2(3 Jo ^ ^ 



(^5) 
(66) 



P 



P 



(67) 



Since ^ea is fixed by Eq. (P^), the second factor on the right-hand-side of Eq. (^) is 
positive (negative) for z_ (-2+). A stable solution corresponds to a negative value of ( |65D 
and (pBD, therefore one has to take the solution q" for Z- and g' for z+. Moreover, since 
for £ < £^TH the right hand side of ( ^4]) is positive then we obtain that C(0) — Pq' and 
C(0) — Pq" are positive and using that 



a > 1 



and imposing that (|66D has to be negative, we obtain: 



1 



Pjp - 1) ^p_2 



> for g = q', q" 



(C{0)-Pq,,y 2 
But this is impossible for z^, i.e. it is not possible to have a consistent stable z^ solution. 
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In this Appendix we give an argument in favor of the relation (3x = da{j3, f)/df in 
Eq. ([611) . Let us start by assuming that it does hold and see that it leads to the equation 
linking x, and C(r) in the Matsubara approach. First of all we write the derivative 
of a with respect to / as a derivative with respect to S. This can be easily done by 
noticing that differentiating / in Eq. (|T^) with respect to £ ai P and T fixed is equivalent 
to differentiating / with respect to S at /?, F, qea and C(r) fixed because / is stationary 
in qea and C(r). As a consequence we find 

da da d£ da , , 

= 57 = a£ a? = a£ ' 

The derivative da/d£ can be easily computed from Eq. (Q). One arrives at 

1 = (C'(0)2-/5V(2;-l)+/3gEAC'(0)(x-2))"' . (69) 

In summary, starting from the TAP approach at fixed £ we obtain and C(r). Assuming 
then that Eq. (|6TD holds we obtain the equation linking ^ba and x in the Matsubara 
approach. The equations for C(r) in the TAP and Matsubara approaches, once ^ea is 
fixed, are identical. Hence we have proven that Eq. (|6lD leads to the Matsubara results 
3§. 



m 



The proof will be complete if we showed the other sense of the implication. 
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